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Detection of Polarimetric Variations Associated with the Shortest 
Time-Scale Variability in S5 0716+714 
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Abstract 

Wc present the result of near-infrared and optical observations of the BL Lac object S5 0716+714 
carried out by the KANATA telescope. S5 0716+714 has both a long term high-amplitude variability and 
a short-term variability within a night. The shortest variability (microvariability) time-scale is important 
for understanding the geometry of jets and magnetic field, because it provides a possible minimum size of 
variation sources. Here, we report the detection of 15-min variability in S5 0716+714, which is one of the 
shortest time-scales in optical and near- infrared variations observed in blazars. The detected microvari- 
ation had an amplitude of 0.061+0.005 mag in V band and a blue color of A{V — J) = —0.025+0.011. 
Furthermore, we successfully detected an unprecedented, short time-scale polarimetric variation which 
correlated with the brightness change. We revealed that the microvariation had a specific polarization 
component. The polarization degree of the variation component was higher than that of the total flux. 
These results suggest that the microvariability originated from a small and local region where the magnetic 
field is aligned. 
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1. Introduction 

Blazars are one of classes of active galactic nuclei 
(AGNs), whose relativistic jets are considered to be di- 
rected along the line of sight (Blandford & Konigl, 1979 
, Antonucci et al., 1993). Blazars show a large and rapid 
variability, and their fiuxes often vary by an order of mag- 
nitude in a time-scale of days. A short time-scale variabil- 
ity has also been reported, that is, so-called intra-day vari- 
ability or microvariability, having a time-scale of < 1 day 
(e.g. Antonucci et al., 1993). 

The observation of the shortest time-scale variability 
gives us important information for understanding the vari- 
ation mechanism, the emitting region and the jet geom- 
etry. An extremely large variability with a time-scale of 
tens of minutes has recently been reported in the opti- 
cal region (e.g. Xie et al., 2004 for eight blazars e.g. 
S5 0716+714). However, Cellone et al., (2007) argued 
that they can be spurious results due to systematic errors 
arising from analysis. The shortest time-scale in blazar 
variability was, in general, believed to be several hours. 

In gamma-ray region, however, a significant variability 
with a time-scale of minutes have recently been discovered 
by Albert et aL, (2007) in Mrk 501 and Aharonian et al., 

(2007) in PKS2155-304. Also in the optical band, a 15- 
min microvariation has been reported by Gupta et al., 

(2008) in Mrk 501. The amplitude of the microvaria- 
tion was not so extremely large, but rather small, with 
0.05 mag in R band. Such a short time-scale indicates that 



the size of the emitting region of the variation is smaller 
than a Schwarzschild radius Rs of the central black hole in 
AGN if we assume typical values for the black hole mass 
and the jet Doppler factor. It is, hence, possible that mi- 
crovariations originate from a small confined region rather 
than a whole emitting area. On the other hand, observa- 
tional features of those shortest time-scale variations have 
not been established due to the lack of detailed obser- 
vations, for example, polarimetric and multi-wavelength 
observations. 

The spectral energy distribution (SED) of blazars 
can be described by synchrotron radiation and inverse- 
Compton scattering. In the optical band, the synchrotron 
radiation is dominant. Variations of polarization can, 
hence, be observed in blazars (e.g. Angel & Stockman, 
1980). The polarimetric observation provides an impor- 
tant clue for the magnetic field and its geometry. However, 
there is still few polarimetric observation of microvariabil- 
ity so far. Cellone et al., (2007) observed hours-scale vari- 
ations in the polarization of AO 0235+164. According 
to their observation, the temporal variation of polariza- 
tion parameters showed no correlation with the variation 
of the total flux. For the shortest time-scale variations 
(~10 min), there was no report studying possible correla- 
tion between the polarization and the total flux. 

S5 0716+714 is a BL Lac object, and one of systems 
in which microvariability was observed. The variability 
has been observed on the time-scale of days and hours 
(e.g. Wagner et al, 1996, Stalin et al., 2006 and Wu et al. 
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2007). Nilsson et al., (2008) estimated its redshift to be 
z = 0.31 ±0.08. 

In order to study the mechanism of the microvariabil- 
ity and to find the shortest time-scale variations, we per- 
formed photopolarimetric observations of S5 0716-1-714 in 
the optical and near-infrared (NIR) bands simultaneously. 
In this paper, we report our detection of 15-min variability 
with polarimetric variations. The 15-min variability has a 
small amplitude, but the light-curve is very well-sampled. 
The paper is arranged as follows: In section 2, we present 
the observation method and analysis. In section 3, we re- 
port the result of the photometry and polarimetry. We 
discuss the emitting region and the geometry of the mag- 
netic field based on the observation in section 4. Finally, 
the conclusion is drawn in section 5. 

2. Observation 

In this paper, we focus on the observation on 20 Oct 
2007. It is part of our monitoring of this blazar in 2007 
using TRISPEC attached to the 1.5-m "KANATA" tele- 
scope at Higashi-Hiroshima Observatory. TRISPEC has 
a CCD and two InSb arrays, producing photopolarimetric 
data in an optical band and two NIR bands simultaneously 
(Watanabe et al., 2005). A unit of the observing sequence 
consisted of successive exposures at four position angles 
of the half-wave plate; 0°, 45°, 22°5, and 67°5. 

The integration times in each exposure were 63, 55 
and 28 sec for V-, J- and ifs-band images, respectively. 
Our one polarimetric data was derived from each set of 
the four exposures. Thus, the time resolution of our 
polarimetry is about 5-min. while that of our photom- 
etry is about 1-min. The sky condition was clear and 
stable during our observation. The full width of half 
maximum of the point spread function was typically 3 
arcsec. All images were bias-subtracted and fiat-fielded, 
and then we performed aperture photometry with the 
IRAF APPHOT package. We adopted differential pho- 
tometry with a comparison star taken in the same frame, 
and the aperture size was about 7 arcsec. Its position 
is R.A.=07^21'°52?3, Dec=-h71°18'17."6 (J2000.0), and 
the magnitudes are 12.48, 11.32 and 10.98 mag in V, 
J, and Ka, bands, respectively (Gonzalez — Perez et al., 
2001; Skrutskie et al, 2006). We confirmed that the 
comparison star kept constant during the observation. 
The dispersion of the magnitudes of the comparison star 
was less than 0.005 mag, calculated with a check star at 
R.A.= 07''21'"54?4, Dec=+71°19'21."3 which is also taken 
in the same frame of S5 0716-1-714, and its magnitudes of 
that are 13.55 and 12.34 in V and J bands. The com- 
parison and the check stars are listed in Ghisellini et al., 
(1997) as star B and C, respectively. 

The instrumental polarization was confirmed to be 
smaller than 0.1 % in F band by observing unpolarized 
standard stars. We applied no correction. The zero point 
of the polarization angle is corrected as the standard sys- 
tem ( measured from north to east ) by observing the 
polarized stars. For NIR bands, the polarimetric accu- 
racy was worse (> 0.4 %) and we used the polarization 
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Fig. 1. Top panel: Light curves of V, J bands and V — J 
color on 20 Oct. 2007. The time bin sizes are 1- and 6-min 
in V and J bands, respectively. Bottom panel: Differential 
magnitudes of the comparison star against the check star in V 
and J bands, indicating that the magnitude of the comparison 
star was constant within 0.005 and 0.010 mag in V and J 
bands, respectively. 

parameter only in V band in this paper. 
3. Result 

3.1. Photometric observation 

Raiteri et al. (2003) reported that the object is at about 
13 mag and 15 mag in V band in its bright and faint states, 
respectively, based on their 8-year light curve from 1994 
to 2002. In Oct 2007, the object was in a bright state 
at 12.6 — 13.5 mag. Figure 1 shows the light curves on 
20 Oct. 2007 when the object was the brightest in the 
period of our monitoring. During 3 hours of our observa- 
tion on 20 Oct., the object brightened by 0.123±0.005 and 
0.101 ±0.010 mag in V and J bands, respectively. In addi- 
tion, there was a bump-like structure around 18.0 (UT) on 
the overall brightening trend. The rising time of the bump 
was about 15-min both in V and J bands. The amplitude 
of the bump was 0.061±0.005 and 0.035±0.010 mag in V 
and J bands, respectively. We did not discuss the light 
curve in A's band, because the error of each photometric 
point is large, about 0.1 mag. 

The time-scale Ar of the bump is estimated as 



1 (Af) 
1 + z dF/dt ' 



(1) 



where 1 + z stands for a cosmological effect (Romero et al., 
2002). The rising and decaying time-scales of the bump 
were 970 and 620 sec in the observer frame. These val- 
ues were estimated after the overall brightening trend was 
subtracted. With z = 0.31±0.08, the rising and decaying 
time-scales are calculated as 740±50 and 480±30 sec in 
the object frame, respectively. 

Cellone et al. (2000, 2007) insisted that the violent 
variations of tens of minutes time-scale may be spurious. 
We now examine the 15-min variation that we detected in 
S5 0716+714 along with the points proposed by Cellone 
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et al. (2000, 2007). First, the spurious variation in the 
hght curve can arise from the variation of the seeing size, 
because of the effect of the host galaxy (Cellone et al., 
2000). However, the host galaxy of S5 0716+714 is so 
faint that its contribution is negligible to our photometry 
(Nilsson et al., 2008). Moreover, we confirmed that there 
is no correlation between the light curve and the tempo- 
ral variation of the seeing size in our case. Second, if the 
comparison star is ^2 mag brighter than the source, the 
variation is severely overestimated (Cellone et al., 2007). 
The magnitude of comparison is 12.48 mag in V band as 
mentioned in section 2, and that of the object was about 
12.6 in the same band. The difference of magnitude in the 
object and the comparison is < 1 mag, thus we did not 
overestimate the variation. In order to confirm the signif- 
icance of the observed bump, we calculate the variability 
parameter C (Romero et al., 2002) with the scale factor F 
(Howell et al., 1988 and Cellone et al, 2007). According 
to Romero et al., (2002). an object can be considered as 
variable (at the 99 % confidence level) if their parameters 
satisfies C/F > 2.576. When the magnitudes of the object, 
comparison and check star were 12.6, 12.48 and 13.55, the 
F was about 0.88, and the scaled variability parameter 
C/F around the bump (t < 17.68 and 18.27 < t in UT) 
was 5.340. Thus the bump was not caused by the effect 
of the magnitudes of the comparison and the check star. 
Therefore, we conclude that the observed 15-min bump is 
a real one. 

The V — J color became bluer with the overall bright- 
ening trend. The color at the end of the observation was 
0.025+0.011 mag bluer than that at the start. Also at 
the 15-min bump maximum, it became bluer, A{V — J) = 
—0.025+0.011, and this behavior is called as "bluer-when- 
brighter" . Such a feature of bluer- when-brighter was also 
observed in the internight variation of one day time-scale 
in S5 0716+714. In a long term trend of tens of days 
time-scale, however, the feature of bluer- when-brighter 
was seldom seen (Wagner et al., 1996 and Ghisellini et al., 
1997). This is the first time that this feature was detected 
in such a short variability (15-min bump) in optical and 
NIR bands. Our observation indicates that microvariabil- 
ity has the same feature as the internight variability. It 
suggests that the internight variability and the microvari- 
ability are caused by the same mechanism, while the long 
term variation has a different one. 

In our monitoring of S5 0716+714 in 2007, there are the 
other two days in which we obtained long time-series data, 
that is, 3 Nov. (3.5 hours) and 13 Nov. (2.5 hours). We in- 
vestigate whether the frequency of microvariability is high 
or low using the scaled variability parameter C/F. For our 
observations, the scaled variability parameters C/F are 
1.22 and 1.04 on 3 and 13 Nov, respectively, in V band. 
Based on this result, we conclude that the microvariability 
only appeared on 20 Oct in our monitoring. 

3.2. Polarimetric observation 

The polarization parameters were variable on 20 Oct, 
too. The maximum polarization degree of the object was 
about 10.5 % and minimum of that was about 7.5 %. 
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Fig. 2. Stokes Q and U parameters in V band during obser- 
vations on 20 Oct. 2007. The Stokes parameters have been 
rcbinned to 4 min. The unit of the Q and U is erg sec~^ cm~^. 

Figure 2 shows the variation on the Q — U plane. Time- 
series of the Stokes parameters, Q and U, are depicted in 
figure 3, with the y-band light curve for comparison. The 
temporal evolution of the Stokes parameters indicates the 
presence of a variation feature on the Q ~ U plane asso- 
ciated with the bump. We calculated differential polar- 
ization vectors from the overall trend. We assumed that 
there were two components in the polarization vector; the 
first component was Phase (0 associated with the overall 
brightening trend and the second one was Pdiff(i) which 
was a deviation from Pbaso(i)- They are written as, 

Pobs(t) = Phase (t) + Pdiff (t) (2) 
Phase (t ) = (Qhase {t) , C/hase (^ ) ) ■ (3) 

In our analysis, Qbasc{t) and C/hase(0 ^^'^ approximated 
as linear functions which are obtained by fitting the data 
outside of the bump {t < 17.68 and 18.27 < t in UT), 
as shown in figure 3. We calculated the differential po- 
larized flux, PFdiff(t) = |Pdiff(t)|, and angle, PAdie{t) = 
0.5arctan[Udis(t) / Qdisit)], which are shown in figure 4. 

If the bump component has no specific polarization vec- 
tor, the PFdiff(<) and PAdm{t) should be mainly con- 
structed by a random noise. In figure 4, however. Pi^diff (0 
is correlated with the bump. The PAdiff(t) dramatically 
changed just before the bump at 17.5 UT, and then, it 
kept constant during the bump. Thus the bump had the 
intrinsic polarization vector whose polarization angle was 
constant. This result is insensitive to the data period 
used for the definition of the overall trend. Although the 
PA-dieit) apparently continues to be constant even after 
the bump, it may be an artificial effect, due to the linear 
fitting to the overall trend. 

Figure 4 indicates another bump-like structure in 
PFdis{t) at 19.5 UT. It can also be associated with a pos- 
sible bump in the light curve, as can be seen from figure 1 
and 3. Thus, such a short bump may frequently appear 
in the brightest state of S5 0716+714. 
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Fig. 3. Top panel: The flux light curve in V band. Middle 
and bottom panel: Q and U curves in the same band. Solid 
lines indicate the Phase which is defined in the text. All 
units of the vertical axis are erg sec~^ cm~^. The flux was 
calculated with 1.98 xlO~^erg sec~^ cm~^ at mag in V 
band (Fukugita ct al., 1995). 
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Fig. 4. Time variation of the difli'erential polarization vec- 
tor from the overall brightening trend. The top flgure shows 
the differential polarized flux P-Fdifl(i) (^rg sec~^ cm~-^). At 
around 18 UT, there is a peak which is associated with the 
bump. The bottom figure shows the differential polarization 
angle PAdiff(t) (degree). 



4. Discussion 

The time-scales in the object frame were 740±50 and 
480±30 sec in the rising and decaying phases of the bump. 
The size of the emitting region, R, can be estimated as 
a light-crossing time, R < cAr, where c is the speed of 
light. With the correction of the Doppler boosting, it 
is now, R < 6cAt. A Schwarzschild radius in the black 
hole mass M is defined as i?g = 2GM/c^, where G is 
the gravitational constant. If the object has typical val- 
ues of the Doppler factor (5wlO and the black hole mass 
M = IO^Mq, the emitting region R is constrained to be 
R < (0.49±0.03)-i?g. In this case, thus, the variation orig- 
inated from a quite confined region whose size is smaller 
than i?g. We now assume that the bump was caused by 
variations of the whole area in the jet at a certain distance 
apart from the black hole. In this situation, we can expect 
that the emitting region is larger than the Schwarzschild 
radius. R > i?g, except for a scenario with an extreme re- 
confinement of the jet (Sokolov et al., 2004). In this case, 
however, the emitting region is smaller than Rg- In order 
to satisfy R> Rs, M should be 

M c^At 



Substituting the decaying time-scale of 480 sec to equa- 
tion (4), we obtain 5 > (20±1) x (M/IO^Mq). The bump 
could, thus, originate from the whole area of the jet if 
S5 0716-1-714 has a relatively small black hole mass or a 
large Doppler factor. 

The whole emitting region in blazars is, however, be- 
lieved to be not so confined in l-Rs- If the emitting 
region is 5 or 10-i?g, the relationship with 6 and Af is 
about 6 > (lOOor2OO)x(Af/lO9M0). In general, the op- 
tical emitting region is believed to be extended in sub-pc 
scale in jets. If the emitting region is located at the sub-pc 
scale from the black hole, the black hole mass (or Doppler 
factor) would be extraordinarily small (or large) , in order 
to explain the bump by the variation of the whole emit- 
ting region. Thus, an alternative idea is required. The 
idea is that the emission region of the microvariability is 
a small and local area in the whole area of the jet. This 
suggestion is consistent with the small amplitude of the 
magnitude at the bump. 

In section 3, we showed that the 15-min bump has 
an intrinsic polarization vector. We estimate the polar- 
ization degree of the bump component. The maximum 
value of the differential polarized flux is (2.7±0.5)xl0~^^ 
erg sec~^ cm~^. Using the total flux at the bump peak, 
(1.02±0.02)xl0~-^^ erg sec~^ cm~^, the polarization de- 
gree of the bump component is calculated to be 27±5%. 
On the other hand, the observed polarization degree at 
the bump is 9.8±0.5%. Therefore, the polarization de- 
gree of the bump component is much larger than that of 
the emission component of the overall brightening trend. 
The bump, hence, presumably originated from the region 
where the magnetic field is more aligned than that in the 
emitting region of the overall brightening trend. 

On the basis of both variability time-scale and the po- 
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larization behavior, we suggest that the bump originated 
not from the whole jet, but from the local region where 
the magnetic field is ahgncd. 

5. Conclusion 

We obtained multicolor photometric and polarimetric 
data of the shortest time-scale variation in S5 0716+714 
for the first time. Our findings are summarized as below; 
first, the object became blue with A{V — J) = —0.025 
during the bump. Second, the bump component has a 
specific polarization vector with the large polarization 
degree of 27%, which is distinct from the overall bright- 
ening trend. It can be suggested that the emitting region 
of the microvariability in optical band is a small and local 
area compared with the whole optical emitting region in 
the jet. 

This work was partly supported by a Grant-in-Aid 
from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan (19740104). 
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